Search for: All records

The water vapor transport associated with latent heat flux (LE) in the planetary boundary layer (PBL) is critical for the atmospheric hydrological cycle, radiation balance, and cloud formation. The spatiotemporal variability of LE and water vapor mixing ratio (r_v) are poorly understood due to the scale‐dependent and nonlinear atmospheric transport responses to land surface heterogeneity. Here, airborne in situ measurements with the wavelet technique are utilized to investigate scale‐dependent relationships among LE, vertical velocity (w) variance (), andr_vvariance () over a heterogeneous surface during the Chequamegon Heterogeneous Ecosystem Energy‐balance Study Enabled by a High‐density Extensive Array of Detectors 2019 (CHEESEHEAD19) field campaign. Our findings reveal distinct scale distributions of LE, , and at 100 m height, with a majority scale range of 120 m–4 km in LE, 32 m–2 km in , and 200 m–8 km in . The scales are classified into three scale ranges, the turbulent scale (8–200 m), large‐eddy scale (200 m–2 km), and mesoscale (2–8 km) to evaluate scale‐resolved LE contributed by and . The large‐eddy scale in PBL contributes over 70% of the monthly mean total LE with equal parts (50%) of contributions from and . The monthly temporal variations mainly come from the first two major contributing classified scales in LE, , and . These results confirm the dominant role of the large‐eddy scale in the PBL in the vertical moisture transport from the surface to the PBL, while the mesoscale is shown to contribute an additional ∼20%. This analysis complements published scale‐dependent LE variations, which lack detailed scale‐dependent vertical velocity and moisture information.

 

Abstract The magnitude of water vapor content within the near-storm inflow can either support or deter the storm’s upscale growth and maintenance. However, the heterogeneity of the moisture field near storms remains poorly understood because the operational observation network lacks detail. This observational study illustrates that near-storm inflow water vapor environments are both significantly heterogeneous and different than the far-inflow storm environment. This study also depicts the importance of temporal variation of water vapor mixing ratio (WVMR) to instability during the peak tornadic seasons in the U.S. Southeast and Great Plains regions during the Verification of the Origins of Rotation in Tornadoes Experiment Southeast 2018 (VSE18) campaign and the Targeted Observation by Radar and UAS of Supercells (TORUS) campaign, respectively. VSE18 results suggest that the surface processes control WVMR variation significantly in lower levels, with the highest WVMR mainly located near the surface in inflows in the southeast region. In contrast, TORUS results show more vertically homogeneous WVMR profiles and rather uniform water vapor distribution variation occurring in deep, moist stratified inflows in the Great Plains region. Temporal water vapor variations within 5-min periods could lead to over 1000 J kg −1 CAPE changes in both VSE18 and TORUS, which represent significant potential buoyancy perturbations for storms to intensify or decay. These temporal water vapor and instability evolutions of moving storms remain difficult to capture via radiosondes and fixed in situ or profiling instrumentation, yet may exert a strong impact on storm evolution. This study suggests that improving observations of the variability of near-storm inflow moisture can accurately refine a potential severe weather threat. Significance Statement It has long been recognized that better observations of the planetary boundary layer (PBL) inflow near convective storms are needed to improve severe weather forecasting. The current operational networks essentially do not provide profile measurements of the PBL, except for the sparsely spaced 12-hourly sounding network. More frequent geostationary satellite observations do not provide adequately high vertical resolution in the PBL. This study uses airborne lidar profiler measurements to examine moisture in the inflow region of convective storms in the Great Plains and the southeastern United States during their respective tornadic seasons. Rapid PBL water vapor variations on a ∼5 min time scale can lead to CAPE perturbations exceeding 1000 J kg −1 , representing significant perturbations that could promote storm intensification or decay. Severe thunderstorms may generate high-impact weather phenomena, such as tornadoes, high winds, hail, and heavy rainfall, which have substantial socioeconomic impacts. Ultimately, by contrasting characteristics of the convective storm inflow in the two regions, this study may lead to a more accurate assessment of severe weather threats. 

The momentum transport by orographic gravity waves (OGWs) plays an important role in driving the large-scale circulation throughout the atmosphere and is subject to parameterization in numerical models. Current parameterization schemes, which were originally developed for coarse-resolution models, commonly assume that unresolved OGWs are hydrostatic. With the increase in the horizontal resolution ofstate-of-the-artnumerical models, unresolved OGWs are of smaller horizontal scale and more influenced by nonhydrostatic effects (NHE), thus challenging use of the hydrostatic assumption. Based on the analytical formulas for nonhydrostatic OGWs derived in our recent study, the orographic gravity wave drag (OGWD) parameterization scheme in the Model for Prediction Across Scales is revised by accounting for NHE. Global simulations with 30-km horizontal resolution are conducted to investigate NHE on the momentum transport of OGWs and their impacts on the large-scale circulation in boreal winter. NHE are evident in regions of complex terrain such as the Tibetan Plateau, Rocky Mountains, southern Andes, and eastern Antarctica. The parameterized surface wave momentum flux can be either reduced or enhanced depending on the relative importance of NHE and model physics–dynamics interactions. The NHE corrections to the OGWD scheme significantly reduce the easterly biases in the polar stratosphere of the Northern Hemisphere, due to both weakened OGWD in the upper troposphere and lower stratosphere and suppressed upward propagation of resolved waves into the stratosphere. However, the revised OGWD scheme only has a weak influence on the large-scale circulation in the Southern Hemisphere during boreal winter.

 

Abstract. High ice water content (HIWC) regions in tropical deep convective clouds, composed of high concentrations of small ice crystals, were not reproduced by Weather Research and Forecasting (WRF) model simulations at 1 km horizontal grid spacing using four different bulk microphysics schemes (i.e., the WRF single‐moment 6‐class microphysics scheme (WSM6), the Morrison scheme and the Predicted Particle Properties (P3) scheme with one- and two-ice options) for conditions encountered during the High Altitude Ice Crystals (HAIC) and HIWC experiment. Instead, overestimates of radar reflectivity and underestimates of ice number concentrations were realized. To explore formation mechanisms for large numbers of small ice crystals in tropical convection, a series of quasi-idealized WRF simulations varying the model resolution, aerosol profile, and representation of secondary ice production (SIP) processes are conducted based on an observed radiosonde released at Cayenne during the HAIC-HIWC field campaign. The P3 two-ice category configuration, which has two “free” ice categories to represent all ice-phase hydrometeors, is used. Regardless of the horizontal grid spacing or aerosol profile used, without including SIP processes the model produces total ice number concentrations about 2 orders of magnitude less than observed at −10 ∘C and about an order of magnitude less than observed at −30 ∘C but slightly overestimates the total ice number concentrations at −45 ∘C. Three simulations including one of three SIP mechanisms separately (i.e., the Hallett–Mossop mechanism, fragmentation during ice–ice collisions, and shattering of freezing droplets) also do not replicate observed HIWCs, with the results of the simulation including shattering of freezing droplets most closely resembling the observations. The simulation including all three SIP processes produces HIWC regions at all temperature levels, remarkably consistent with the observations in terms of ice number concentrations and radar reflectivity, which is not replicated using the original P3 two-ice category configuration. This simulation shows that primary ice production plays a key role in generating HIWC regions at temperatures <-40 ∘C, shattering of freezing droplets dominates ice particle production in HIWC regions at temperatures between −15 and 0 ∘C during the early stage of convection, and fragmentation during ice–ice collisions dominates at temperatures between −15 and 0 ∘C during the later stage of convection and at temperatures between −40 and −20 ∘C over the whole convection period. This study confirms the dominant role of SIP processes in the formation of numerous small crystals in HIWC regions. 

Planetary boundary layer (PBL) schemes parameterize unresolved turbulent mixing within the PBL and free troposphere (FT). Previous studies reported that precipitation simulation over the Amazon in South America is quite sensitive to PBL schemes and the exact relationship between the turbulent mixing and precipitation processes is, however, not disentangled. In this study, regional climate simulations over the Amazon in January–February 2019 are examined at process level to understand the precipitation sensitivity to PBL scheme. The focus is on two PBL schemes, the Yonsei University (YSU) scheme, and the asymmetric convective model v2 (ACM2) scheme, which show the largest difference in the simulated precipitation. During daytime, while the FT clouds simulated by YSU dissipate, clouds simulated by ACM2 maintain because of enhanced moisture supply due to the enhanced vertical moisture relay transport process: (a) vertical mixing within PBL transports surface moisture to the PBL top, and (b) FT mixing feeds the moisture into the FT cloud deck. Due to the thick cloud deck over Amazon simulated by ACM2, surface radiative heating is reduced and consequently the convective available potential energy is reduced. As a result, precipitation is weaker from ACM2. Two key parameters dictating the vertical mixing are identified,p, an exponent determining boundary layer mixing andλ,a scale dictating FT mixing. Sensitivity simulations with alteredp,λ, and other treatments within YSU and ACM2 confirm the precipitation sensitivity. The FT mixing in the presence of clouds appears most critical to explain the sensitivity between YSU and ACM2.

 

Abstract Observation system simulation experiments are used to evaluate different dual-Doppler analysis (DDA) methods for retrieving vertical velocity w at grid spacings on the order of 100 m within a simulated tornadic supercell. Variational approaches with and without a vertical vorticity equation constraint are tested, along with a typical (traditional) method involving vertical integration of the mass conservation equation. The analyses employ emulated radar data from dual-Doppler placements 15, 30, and 45 km east of the mesocyclone, with volume scan intervals ranging from 10 to 150 s. The effect of near-surface data loss is examined by denying observations below 1 km in some of the analyses. At the longer radar ranges and when no data denial is imposed, the “traditional” method produces results similar to those of the variational method and is much less expensive to implement. However, at close range and/or with data denial, the variational method is much more accurate, confirming results from previous studies. The vorticity constraint shows the potential to improve the variational analysis substantially, reducing errors in the w retrieval by up to 30% for rapid-scan observations (≤30 s) at close range when the local vorticity tendency is estimated using spatially variable advection correction. However, the vorticity constraint also degrades the analysis for longer scan intervals, and the impact diminishes with increased range. Furthermore, analyses using 30-s data also frequently outperform analyses using 10-s data, suggesting a limit to the benefit of increasing the radar scan rate for variational DDA employing the vorticity constraint. 

A suite of six idealized supercell simulations is performed in which the surface drag coefficient C_dis varied over a range of values from 0 to 0.05 to represent a variety of water and land surfaces. The experiments employ a new technique for enforcing a three-force balance among the pressure gradient, Coriolis, and frictional forces so that the environmental wind profile can remain unchanged throughout the simulation. The initial low-level mesocyclone lowers toward the ground, intensifies, and produces a tornado in all experiments with C_d≥ 0.002, with the intensification occurring earlier for larger C_d. In the experiment with C_d= 0, the low-level mesocyclone remains comparatively weak throughout the simulation and does not produce a tornado. Vertical cross sections through the simulated tornadoes reveal an axial downdraft that reaches the ground only in experiments with smaller C_d, as well as stronger corner flow in experiments with larger C_d. Material circuits are initialized enclosing the low-level mesocyclone in each experiment and traced backward in time. Circulation budgets for these circuits implicate surface drag acting in the inflow sector of the supercell as having generated important positive circulation, and its relative contribution increases with C_d. However, the circulation generation is similar in magnitude for the experiments with C_d= 0.02 and 0.05, and the tornado in the latter experiment is weaker. This suggests the possible existence of an optimal range of C_dvalues for promoting intense tornadoes within our experimental configuration.

 

Ensemble Kalman filter (EnKF) analyses of the storms associated with the 8 May 2017 Colorado severe hail event using either the Milbrandt and Yau (MY) or the NSSL double-moment bulk microphysics scheme in the forecast model are evaluated. With each scheme, two experiments are conducted in which the reflectivity ( Z) observations update in addition to dynamic and thermodynamic variables: 1) only the hydrometeor mixing ratios or 2) all microphysical variables. With fewer microphysical variables directly constrained by the Z observations, only updating hydrometeor mixing ratios causes the forecast error covariance structure to become unreliable, and results in larger errors in the analysis. Experiments that update all microphysical variables produce analyses with the lowest Z root-mean-square innovations; however, comparing the estimated hail size against hydrometeor classification algorithm output suggests that further constraint from observations is needed to more accurately estimate surface hail size. Ensemble correlation analyses are performed to determine the impact of hail growth assumptions in the MY and NSSL schemes on the forecast error covariance between microphysical and thermodynamic variables. In the MY scheme, Z is negatively correlated with updraft intensity because the strong updrafts produce abundant small hail aloft. The NSSL scheme predicts the growth of large hail aloft; consequently, Z is positively correlated with storm updraft intensity and hail state variables. Hail production processes are also shown to alter the background error covariance for liquid and frozen hydrometeor species. Results in this study suggest that EnKF analyses are sensitive to the choice of MP scheme (e.g., the treatment of hail growth processes).

 

Since the 1950s, a countergradient flux term has been added to some K-profile-based first-order PBL schemes, allowing them to simulate the slightly statically stable upper part of the convective boundary layer (CBL) observed in a limited number of aircraft soundings. There is, however, substantial uncertainty in inferring detailed CBL structure, particularly the level of neutral stability (zn), from such a limited number of soundings. In this study, composite profiles of potential temperature are derived from multiyear early afternoon radiosonde data over Beijing, China. The CBLs become slightly stable above zn ~ 0.31–0.33zi, where zi is the CBL depth. These composite profiles are used to evaluate two K-profile PBL schemes, the Yonsei University (YSU) and Shin–Hong (SH) schemes, and to optimize the latter through parameter calibration. In one-dimensional simulations using the WRF Model, YSU simulates a stable CBL above zn ~ 0.24zi, while default SH simulates a thick superadiabatic lower CBL with zn ~ 0.45zi. Experiments with the analytic solution of a K-profile PBL model show that adjusting the countergradient flux profile leads to significant changes in the thermal structure of CBL, informing the calibration of SH. The SH scheme replaces the countergradient heat flux term in its predecessor YSU scheme with a three-layer nonlocal heating profile, with fnl specifying the peak value and z*SL specifying the height of this peak value. Increasing fnl to 1.1 lowers zn, but to too low a value, while simultaneously increasing z*SL to 0.4 leads to a more appropriate zn ~ 0.36zi. The calibrated SH scheme performs better than YSU and default SH for real CBLs.